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INTRODUCTION 


There  are  many  aspects  to  the  design  of  rifling  in  guns.  Here,  we  are  concerned  with 
only  one -specification  of  the  rifling  angle  or  slope  to  produce  a  desired  projectile  torque.  A 
rigid  gun  and  a  pre-engraved  projectile  are  assumed.  No  frictional,  thermal,  or  plastic  effects  are 
considered. 

Constant  twist  rifling  has  been  in  use  ever  since  the  earliest  rifled  carmon,  and  in  early 
cannon  which  used  pre-engraved  projectiles,  constant  twist  rifling  was  mandatory.  With  the 
advent  of  the  soft  metal  rotating  band,  however,  it  became  possible  to  vary  the  rifling  angle  for 
various  reasons. 

Progressive  twist  rifling,  for  instance,  can  be  used  to  reduce  projectile  torque  near  the 
origin  of  rifling  where  the  base  pressure  is  highest,  while  constant  twist  rifling  produces  a  torque 
which  is  proportional  to  the  base  pressure.  Progressive  twist  can  be  particularly  important  after 
wear  has  taken  place  near  the  origin  of  rifling.  Progressive  twist  can  still  engrave  the  rotating 
band  reliably,  while  the  velocity  gained  by  the  projectile  over  the  worn  section  of  a  constant  twist 
gun  can  ruin  the  rotating  band  as  it  hits  the  unworn  area. 

Unfortunately,  most  efforts  at  variable  twist  rifling  design  have  used  analytical  forms  for 
the  rifling  curve  and  then  computed  the  resulting  torque  in  retrospect.  This  trial  and  error 
method  may  have  been  convenient  before  the  advent  of  computers,  but  it  is  now  an  unnecessarily 
restrictive  design  technique.  For  instance,  although  traditional  progressive  twist  rifling  does 
decrease  torque  near  the  origin  of  rifling,  it  raises  the  torque  considerably  near  the  muzzle, 
whereas  constant  twist  rifling  has  a  considerably  lower  torque  near  the  muzzle  due  to  the  lower 
base  pressure  there. 

It  would  seem  that  we  cannot  get  everything  that  we  want,  namely  lower  torques  near  the 
origin  of  rifling  and  no  torque  at  all  as  the  projectile  exits  the  muzzle.  Fortunately,  what  seems 
to  be  the  case  is  not  the  case  at  all.  In  fact,  if  we  are  prepared  to  spend  a  few  seconds  of 
computer  time  doing  some  simple  numerical  integrations,  we  can  have  virtually  any  shape  torque 
curve  we  want.  Another  comment  is  perhaps  in  order  for  the  uninitiated;  although  the  purpose 
of  rifling  and  the  torque  produced  thereby  is  to  spin  stabilize  the  projectile  in  flight  for  far 
greater  accuracy,  the  shape  of  the  torque  curve  has  absolutely  nothing  whatsoever  to  do  with  the 
amount  of  spin  imparted  to  the  projectile  in  flight.  The  only  determinant  of  this  spin  is  the  slope 
of  the  rifling  curve  at  the  muzzle.  As  long  as  this  slope  at  this  single  point  is  held  constant,  the 
exit  spin  will  remain  constant  and  we  are  free  to  pick  virtually  any  shape  torque  curve  we  want. 
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RIFLING  SYNTHESIS 


Let  P{x)  be  the  pressure  travel  curve,  A  be  the  bore  area,  r  be  the  bore  radius,  and  m  be 
the  projectile  mass.  The  movement  of  the  projectile  through  the  rifled  bore  is  described  by  the 
following  equations: 

PA-nSin4>=m  — 
dt 

L=nrCos<t>=i— 

dt^ 

where  n  is  the  collective  normal  force  on  the  lands,  0  is  the  rifling  angle,  6  is  the  amount  of 
rotation  of  the  projectile,  and  i  is  the  projectile  moment  of  inertia.  L  is  the  torque  applied  to 
the  projectile  and  v  is  its  velocity. 

Let  f{x)  be  the  rifling  curve  which  is  related  to  the  amount  of  rotation  by 

m=rd(x) 


First, 


and 


^=v— 
dt  dx 


=-f(x) 


dt^ 


V  d 
r  dx 


(vf'ix)) 


dv  _^  dv  _1  d  2 
dt  dx  2  dx 


Hence,  the  equations  of  motion  become 

AP{x)  -  nSinip—m—v^ 
2  clx 

L=nrCos<p=—-^(yf^(x)) 
r  dx 


Solving  for  n  in  the  second  equation  gives  us 


iv  d 

r^Cos<pdx 


iyf'ix)) 
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Inserting  this  in  the  first  equation  yields 


APix)-^Tan<t>-^(yf'{x))= 

ax  ^  dx 


but 


Tan<l>=f'(x) 


Hence,  we  have 


AP(x)  =  ^J-A(yf'(x)f 

^  dx 


or 


m  ax 


*—Ayf'(x))^ 

mr^  j 


Integrating  this,  we  have 


where  R  is  the  projectile  radius  of  gyration.  Hence,  projectile  velocity  is  given  by 


v(x)  = 


Let  /  be  the  length  of  travel  (of  the  rotating  band).  P(x)  is  the  projectile  base  pressure  when  the 
rotating  band  is  at  x.  A  good  typical  value  for 

fy(l) 

r 

is  0.1.  We  therefore  see  that  the  muzzle  velocity  of  a  rifled  gun  is  about  one-half  of  one  percent 
less  than  the  smoothbore  muzzle  velocity.  Indeed  we  can  see  as  Finlayson  points  out  that  the 
projectile  velocity  is  roughly  invariant  with  respect  to  virtually  all  forms  of  rifling  (ref  1). 
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Now,  define  the  torque  function  to  be  the  product  of  some  constant  and  some  torque 
shape  function: 


Ux)  =  cX{x) 


Recalling  that 


A(y(x)f'(x))  = 


rUx) 

iv(x) 


We  have 


4-(^(x)f'(x))  = 
ax 


cr  X(x) 
i  v(x) 


and  integrating  this  equation,  we  get 


Now  let 


f(x)  = 


i  v(x)Jo  v(u) 


I{X)- 


£Piu)du 


and 


Hence,  v  is  given  by 


vix)  = 


2A  I(x) 


\  m  U(x) 


and  therefore 


f/(x)  =  nMu)U(u)^ 

■'  2Ai  I(x)Jo  i(u) 
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Letting 


JH  I(u) 


We  have 

f/M  - 

^  2AR^  I(x) 

And  hence,  when/’(/)  has  been  specified,  c  is  defined  as 

^_2ARYmi) 

rUmi) 


Therefore 

w.,.f(im)U(x)J(x) 
^  U(l)J(l)I(x) 


Rewriting  this  last  equation  in  the  form 

f(x)=f(l)^g(x) 


where 


we  have 


g(x)  = 


imx) 

mix) 


which  we  can  solve  for/’(x),  getting 

f'(x)  =  - 


f'iMx) 


m-fWn-gixf) 

r 


5 


Now,  for  the  sake  of  numerical  accuracy  in  computing  the  integral  J(x)  (especially  for 
small  values  of  a:),  it  is  necessary  to  perform  an  integration  by  parts  on  J.  Recall  that  J  is  given 
by 


■'0  I(u) 


We  must  have 


X(0)=0 


for  this  integral  to  exist,  since  /  is  proportional  tox  for  small  values  of  a:. 
Since 


l(x)^=  rP{u)du 
Jo 


we  have 


Hence, 


■'W  .rmmu) 

Jo  P(u)  Jo 


where 


2X(x)U(x) 

P(x) 


Note  that 


K(0)  = 


2X'(0)U(0) 

P'(0) 


by  THospital’s  rule. 

Using  integration  by  parts  on  J,  we  have 

J(x)  =  K(x)I(x)  -  K(0)I(0)  -S(x)= K(x)I(x)  -  S(x) 


where 


Six)=fJl(u)dK(u) 


The  Stieltjes  integral  S  should  be  computed  directly  by  trapezoidal  rule  because  K  may  not  be 
differentiable.  Indeed,  K  need  not  even  be  continuous. 

Now,  in  order  to  compute  g(0),  it  is  sufficient  to  note  the  asymptotic  behavior  of  A,  P,  /, 
and  J  for  small  values  of  x: 


A(x)-xAl[0) 

P(x)-xP''(0) 

I(xf=J^P(u)du  ~  ^x^P'm 


J(x)-  xX'(0)U(0) 


N 


2 

P'(0) 


We  therefore  have 


m  .  2A!mM=K(0) 

I(x)  P'(0) 


and 


We  now  summarize  in  one  place  the  equations  applicable  to  rifling  synthesis. 

P{0)  =  0  X{0)=0 


lix)=^fJP(u)du 


U(x)= 

K(x) 

P'(0) 

S(x)=p(u)dK(u) 

J(x)=K(x)I(x)-S(x) 


_  2X(x)U(x) 
P(x) 
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g(x)- 


mm 

jmx) 


g(o)=^m 

j(i) 

f(x)=f(l)^8(x) 


f'ix)  = 


f(l)g(x) 


N 


(l-g(xf) 


m=p'{u)du 


The  equations  defining /(x)  do  so  implicitly  because /(jc)  enters  the  right-hand  sides 
through  U.  Functional  iteration  is  used  to  obtain /'(*).  Convergence  is  extremely  rapid,  taking 
only  a  few  iterations.  In  fact,  if  we  take /’(*)=/’(/)  initially,  the  very  first  iteration  gives  a  result 
which  is  in  error  by  only  about  one  quarter  of  a  percent  for  most  reasonable  A’s.  The  reason  for 
this  is  that  U  is  only  very  weakly  dependent  onf{x),  since  the  square  of  the  latter  is  small 
relative  to  unity.  Hence,  U  is  roughly  a  constant  near  unity  for  all  practical  rifling  configurations. 

Also  note  from  the  equations  that  for  the  first  iteration,  17  is  a  constant  which  cancels 
out.  Hence  the  rifling  slope  function  is  determinable  within  about  a  quarter  of  a  percent  error 
without  knowledge  of  either  the  projectile  radius  or  radius  of  gyration.  The  absolute  values  of 
pressure  are  also  unnecessary.  Only  the  shapes  of  the  pressure  travel  function  and  the  desired 
torque  function  are  necessary  to  determine  a  good  approximation  to  the  rifling  configuration  for 
any  size  gun  with  any  muzzle  twist  and  any  projectile  regardless  of  its  mass  or  moment  of  inertia! 


FIGURES  OF  MERIT 

Let  us  take 


X(x)=kP{x) 


For  the  first  (and  last)  iteration,  set 


f'(x)^f'(l) 


Hence, 

U(x)= const  =U 
K(x)=2kU=const 

S(x)=£‘ldK=0 

J(x)  =  2kUI(x) 

I(l)2kUI(x)^j 
*  2kUimx) 

Therefore 


(immediate  convergence  to  constant  twist  rifling). 


Now,  recall  that 


Torque  =L=cX 


where 

2ARY(1)I(1) 

rU(im 

Denote  the  max  norm  by  /  S.  The  maximum  torque  force  will  then  be 

M=sM 

r  r 


and  the  maximum  pressure  force  is  IP  HA.. 

We  define  /i  to  be  the  ratio  of  maximum  torque  force  to  maximum  pressure  force: 

_  cjAi  _  Ur 
iPiAr  ipmim 
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For  the  constant  twist  case  (AocP),  we  have 


lP/U(l)2kU(l)I(l) 


We  can  also  define  another  figure  of  merit  (a)  as  the  ratio  of  maximum  torque  to  maximum 
constant  twist  torque: 


_2IA/I(l)U(l) 

Pc  IPIM 

Note  that  ct=1  in  the  constant  twist  case  and  that  a  is  virtually  independent  ofR,  r,  and /’(/). 

We  can  now  compare  different  torque  shape  functions  by  comparing  their  fi’s  and  a’s. 

We  would  presumably  like  these  figures  of  merit  to  be  acceptably  small  (up  to  a  point). 
Now,  the  smallest  values  for  p  and  a  one  can  typically  get  (depending  on  the  pressure  travel 
curve  shape)  are  about  0.03  and  0.45,  respectively.  We  get  these  rather  small  values  by  letting  A 
be  constant  except  in  a  small  neighborhood  of  ac=0  where  A'  must  be  quite  large.  Unfortunately, 
the  byproduct  of  this  A  is  an/  having /(O)  considerably  larger  than/'(/),  an  obviously  impractical 
situation.  However,  if  we  let  A  be  proportional  to  P  up  to  maximum  pressure  and  constant 
everywhere  else  except  for  the  last  five  percent  of  travel  where  we  drop  A  linearly  down  to  zero, 
we  get  a  ^  of  about  0.04  and  a  <r  of  about  0.6  with  an/'(0)  of  about  0.6  f{l).  The  rifling  is 
therefore  progressive,  but  not  radically  so,  and  it  substantially  reduces  maximum  torque  on  the 
projectile  while  giving  it  a  torqueless  exit. 


PRESSURE  TRAVEL  FUNCTIONS 

Ideally,  one  might  want  to  work  with  the  best  available  pressure  travel  data  for  a  given 
round;  however,  the  need  will  generally  exist  to  fire  more  than  one  type  of  round  or  charge  in 
the  gun.  If  each  round  type  displayed  the  same  shape  pressure  travel  function,  the  rifling 
configuration  could  be  determined  uniquely.  In  general,  however,  the  pressure  travel  function 
will  vary  in  shape  for  rounds  of  different  types  and  it  will  be  necessary  to  use  some  representative 
or  average  pressure  travel  shape  function.  This  being  the  case,  it  makes  sense  to  approximate 
the  pressure  travel  data  with  some  smooth  closed  form  function. 

If  the  reader  has  some  favorite  numerical  method  for  approximating  pressure  travel  data, 
this  section  can  be  completely  ignored.  In  any  case,  no  claim  of  superiority  is  made  for  the 
simple  approximation  offered  here;  it  at  least  suffices  for  testing  the  software  associated  with  the 
problem  at  hand. 
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A  simple  approximation  one  might  use  is 


P{x)=axe 

To  determine  the  three  parameters  a,  b,  and  c,  one  might  enforce  the  three  conditions 

P(Xi)=Pi 

P'(xj)=0 

P(X2)=P2 

This  would  allow  us  to  put  the  peak  pressure  (pj)  where  we  wanted  it  (xj)  and  specify  the 
pressure  (pj)  at  the  muzzle  (x^). 

This  representation  seems  insufficiently  flexible,  however,  so  the  following  piecewise 
exponential  representation  is  suggested: 

Pj(x)  if  0  £  X  ^  X2 
f(^>- if  X,  .  . 

where 

i-T 

Pj(x)  =  aixe 

1—1 

P2(x)  =  a2+b2e  '  ' 

and  the  six  parameters  aj,bj,Cj,a2,b2<  and  Cj  are  determined  by  the  conditions 

Pi(h)=Pi 

P^(xj)=0 

P 1(^2^  ~P2 
P2(X2)=P2 

Pj{X2)  =  P2iX2) 

P2(^^=P3 
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Determining  the  parameters  in  P]{x): 


i-T 

Pj(x)=ajXe 


Piix)=aj(l 


Pi(xt)=0 


bj=XjCj\ 


I  _\f 


and 


PMi)=Pi 


We  therefore  have 


in  which  we  must  now  determine  Cj. 

We  have  the  order  relationships 

^  ^2  ^  Pl>  P2^  Pi 


Therefore 


P l{^‘^~P2 


^2Pl 
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or 


Letting 


Therefore,  since /(0)=0,/'(0)  <  0,f{x)  >  0  and/(oo)=oo,  a  positive  zero  of/is  guaranteed. 
Also,  there  must  be  an  m  such  that/’(m)=0.  Hence, 


a"lna+0  =  O ,  m  = 

In  a 
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Expanding /in  a  Taylor  series  around  m,  we  have 


Hence, 


f(z)=0 


-2f(m) 


N  f'U 


This  is  our  first  estimate  of  z  (or  Cj)  which  we  will  subsequently  use  in  Newton  iteration  to  obtain 
z  exactly. 


Now, 


m=- 


Ink 

Ina 


Ina 

f"(m)  =  a”(ln  a)^=e”‘"‘'(ln  af=k(ln  af 
f''(m)  (Inaf 


Therefore  the  first  estimate  of  z  reduces  to 


:.-^(lnk^p(lnk^-l)) 
Ina  »  * 


And  the  Newton  iteration  for  z  is 

f(z)  a^+Pz-1  _e^(zlna-l)  +  l 

f'(z)  d^lna+P  a^lna+P 
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Now,  determining  the  parameters  of  P2{x)\ 


So, 


and 


P/x)=a^^b,e  \  I 


p'(x,)=Pj'(x2)=p'  -  b2=-C2P2 


P 2(^2)=  Pi  "*  02=P2~^2=P2*^2P2 


Therefore,  we  have 


P2M=P2*^2P2 


Now, 


6=Xj-X2>0  .  P2(Xj)=Pj  -• 

CiP^i^-^  ^^*P2-P3=^ 


which  must  be  solved  for  c^. 
Define  f  by 


f(x)=xp2 


^Pz-Pi 


Clearly, 


fic2)=0  .  f(0)=P2-Pi  >  0 
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and  for  large  values  of  x, 


X 


Therefore, 


Now, 


f(«>)=dp2+P2-P} 


and  since. 


we  have 


lim  e_ 
x^O 


=0 


f'(0)=P2  <  0 


Also,  for  large  values  of  x. 


6 


e  * 


-  f'M=0 


Furthermore, 


f^x)=  - 


^  0 


Therefore,  f  is  monotone  increasing  and  can  never  be  positive.  We  can  therefore  finally 
conclude  that  /  is  monotone  decreasing  and  is  guaranteed  to  have  exactly  one  positive  zero 
provided 


f{«)=  6p2-<-p2-P3  <  0 
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The  Newton  iteration  for/(2)=0  is 


Taking  2=0  initially,  we  have 


z  « 


m  _Pi-P2 

f(0)~  P2 


as  our  first  guess  at  2. 

Computing 

zf(z)  -f(z)  =P3  -P2  -  ^P2^  * 


we  have  the  Newton  iteration  for  the  correct  2: 


z- 


Pi~P2  ^Pz^ 


Pz 


) 


Summarizing  the  results  of  this  section: 

PAx)  if  0  ^  X  ^  X2 

^X2<X.X, 

g(x)  =  l-(^p 

h 


Pj(x)  =  -g(x)e 


17 


H  I  ft*  M  Ift* 


^2 

a=— 


k= 


K^zPi) 


_  -P 


Ina 


a‘(c, Ina -1)^1 
a‘‘lna*p 


(»■?) 


P2(^)"P2*<^2P2 
P2=Pl(X2> 

S  =  Xj-X2 

need  P2  <  El-El 


C2- 


JrP2 


P2 


P3-P2-^P2^ 


P2 
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RIFLING  ANALYSIS 


For  the  sake  of  completeness,  we  briefly  derive  the  equations  to  compute  torque  and 
pressure  in  terms  of  rifling  twist.  Recalling  that 

L  =!lA(yAx)) - 

r  dx  2rf'(x)  dx 

i  d 


2r/(jc)  ^ 


(v/W 


and  further  recalling  that 


v(Jt)=, 


lA  I(x) 
\  m  U(x) 


we  have 


Ux)= 


i  d 

2A/(*)^/(jc)* 

2rfXx)  dx 

»»  Uixf 

AR^  d 


f'(xff^"P(u)du 


rf'ix)  dx 


1  + 


^Rf!^ 

K  r  j 


For  constant  twist  rifling,  where  f  (A:)=constant=/t,  we  have 

JckrR^ 


Ux)= 


.(MJ 


If  one  wished  to  write  P  in  terms  of  L,  one  could  begin  with 

2A  dx 


and  since 


■^(v/(x))=-^ 

dx  IV 
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which  gives  us 


and  hence  that 


dx  i 

=  j-(yf'(.x)f 
dx 

vyW=^  [^lXu)fXu)du 

I  Jo 


m_±t!M.(‘uu-^,u)du 

2Adx  i  ^ 


r  d\(  1  ,R^ 

AR^dx^/Xxf 


^^Uu)f'(u)du 

0 
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